Correcting detection error in quantum computation and state engineering 

through data processing 
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Quantum error correction in general is experimentally challenging as it requires significant expan- 
sion of the size of quantum circuits and accurate performance of quantum gates to fulfill the error 
threshold requirement. Here we propose a method much simpler for experimental implementation 
to correct arbitrary detection errors. The method is based on processing of data from repetitive 
experiments and can correct detection error of any magnitude, as long as the error magnitude is 
calibrated. The method is illustrated with its application to detection of multipartite entanglement 
from quantum state engineering. 
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A celebrated achievement in quantum information sci- 
ence is establishment of the error threshold theorem, 
which assures that any experimental error below a certain 
threshold value can be corrected in quantum computa- 
tion or communication [l], Q • This theorem is important 
as it indicates there is no fundamental obstacle to real- 
ization of quantum computation with imperfect experi- 
mental devices, as long as the imperfection is small. The 
experimental realization of fault-tolerant quantum error 
correction, unfortunately, is still challenging. When all 
the experiment devices are subject to errors as it is the 
case for real applications, realization of fault-tolerance 
requires significant expansion of the size of quantum cir- 
cuits for complicated encoding and accurate performance 
of quantum gates to fulfill the error threshold require- 
ment, which is still beyond the capability of current ex- 
perimental technology. 

In this paper, we propose an alternative method to 
correct a special but practically important type of error 
- the detection error — in quantum information pro- 
cessing. The detection error is a common source of noise 
in many quantum information experiments. In particu- 
lar, it is a significant obstacle to observation of multi- 
partite entanglement in quantum state engineering Q. 
We show here that this type of error can be corrected 
at any magnitude as long as the error magnitude has 
been calibrated (for instance, through prior test exper- 
iments). The detection error distorts the experimental 
data by a transformation that depends on the magni- 
tudes of various error possibilities. When the relevant 
error magnitudes have been calibrated for the detectors 
by the prior test experiments, the form of the distortion 
transformation induced by the detection error is known, 
and then we can find a way to inverse this transformation 
to reconstruct the original signal. In this way, we can use 
imperfect detectors to simulate perfect detectors as long 
as their imperfection has been calibrated. The proposed 
method is straightforward for experimental implementa- 
tion as it is based on data processing and removes the 
difficulty associated with fault-tolerant quantum encod- 



ing. To correct the detection error, we only require to 
repeat the same experiments by some additional rounds 
to have small statistical error for the inverse transforma- 
tion. To illustrate its applications, the method is used to 
significantly improve the detection of multi-qubit entan- 
glement and spin squeezing. In many cases, the signal of 
multipartite entanglement only becomes visible after the 
proposed correction of the detection error, in particular 
when the number of qubits is large. 

Any measurements in quantum information can be re- 
duced to population measurements in certain bases (in- 
cluding possibly several complementary bases). If we 
want to measure properties associated with a state p 
(generally mixed) of n qubits, in each chosen measure- 
ment basis, there are 2™ possible measurement outcomes. 
By measurements we determine the probability /$ associ- 
ated with each outcome i (i = 1, 2, • • • , 2"). For instance, 
if we repeat the same experiment N times and get the 
ith outcome Ni times, we estimate the probability fi by 
fi = N.-JN and its standard deviation (the error bar) by 
Afi = \J fi(l — fi)/N using the binomial distribution. If 
the detectors are perfect, the measured probabilities fi 
just give the distribution g$ = of the state p in 

the measurement basis {]«)}. In reality, however, the de- 
tectors always have errors, which distort the distribution 
<7i, making the measured distribution fi significantly dif- 
ferent from gi. The purpose of this paper is to show how 
to reconstruct the real distribution gi from the measured 
distorted probabilities fi. 

We first consider the case where the measurements 
haven individual addressing, and each qubit is measured 
by an independent detector. For detection on a qubit, 
the most general error model is characterized by a 2 x 2 
matrix 



D 



1 - Po Pi 
Po 1 - Pi 



(1) 



where po (pi) denotes respectively the error probability 
that the detector gives outcome 1 (0) for the input signal 
of (1). For simplicity of notation, we assume the error 
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matrix D has the same form for detection of each qubit 
(it is straightforward to generalize the formalism to the 
case where the error rates p and pi in the D matrix 
are qubit-dependent). Furthermore, we assume po and 
Pi have been well calibrated by a prior test experiment. 
For instance, we may input a known state to the detector 
and can calibrate po and p\ easily from the measurement 
data. 

For n qubits, the error model for the detection is then 
characterized by a 2™ x 2" matrix M = [Mji], with the 
element Mji corresponding to the probability of record- 
ing the outcome j with the input signal i. Assume the 
detection error rates on different qubits are independent 
to each other and the binary string i has no zeros and 
m = n — n ones. If we need a flips from to 1 and 
j3 flips from 1 to to change the string from i to j, the 
matrix element Mji is given by 



M 



(l-po) n °- a (l-Po) 



ni-0 



Po pi- 



rn 



The measured probabilities fj are connected with the 
real distribution gi through the distortion transformation 
fj = J2i=i Mjigi.To reconstruct the real signal gi from 
the measured distribution fj , in principle we only need to 
inverse the matrix M. However, as M is a huge 2™ x 2" 
matrix, it is not clear how to inverse this matrix (it is 
even a question whether the inverse exists). 

Our key observation is that the matrix M . with the 
elements given by Eq. (2), has a simple tensor product 
structure. It is straightforward to show by mathematical 
induction that 



M = 0D fe , 



(3) 



fe=i 



where all the are identical and given by D in Eq. (1). 
Therefore, the inverse can be easily done in an analytic 
form with 



M~ 



)Dl 



k=i 



k=i 



1 



Po Pi 



(4) 



where the parameters p' and p\ are given by 



Po = Po/(Po+Pi - 1), 
Pi = Pi/(po+Pi - !)• 



(5) 



Note that with the substitution in Eq. (5), M _1 and M 
have the same form except that p' and p[ can not be 
interpreted as error rates any more since in general they 
are not in the range [0, 1]. The formula also shows that 
the inverse transformation M _1 always exists except for 
the special case with po + Pi = 1 • 

In some experimental systems we do not have the abil- 
ity to resolve individual qubits. Instead, we perform col- 
lective measurements on n qubits by detecting how many 
qubits (denoted by j, j = 0, 1, ■ • ■ , n) are in the state |1) 



in a chosen detection basis (this is equivalent to measure- 
ment of the collective spin operator along a certain direc- 
tion) . In this case, the detection only has n + 1 outcomes 
for an n-qubit system. For collective measurements on 
n qubits, the detection error matrix is represented by an 
(n + 1) x (n + 1) matrix L = [L^]. The matrix element 
Lij corresponds to the probability to register outcome i 
when j qubits are in the |1) state. If the detection error 
matrix for an individual qubit is still given by D in Eq. 
(1), we can directly calculate Lij from D: from signal j 
to i, if nio qubits flip from to 1 and noi qubits flip from 
1 to 0, with the constraints < noi < j, < n w < n — j 
and n i - n 10 = j - i, is given by 

L ij = B (3,Pi,noi)B{n - j,p ,n 10 ) 

"01, "10 

= ^ J B{j,l-p 1 ,q)B{n-j,po,i-q), (6) 
<? 

where B(n,p, k) = (™)p fe (l — p) n ~ k and we have let q 
i — mo, and hence q satisfies the constraint max{0, i+j — 
n) < q < min{i,j}. As the dimension of the L matrix 
depends linearly on the qubit number n, it is typically not 
difficult to numerically calculate its inverse matrix L^ 1 
if n is not very large. As will be shown in the appendix, 
there is also a simple analytic formula for L^ 1 = [Lr. j : 
if we denote the dependence of in Eq. (6) on po,Pi 
as = (po,pi), we have 



(7) 



where p' ,p'i are given by the simple substitution in Eq. 
(5). With the inverse matrix L^ 1 , the real signal g t can 
be similarly reconstructed from the measured data fj as 

9i = ^2j fj- 

The above formulation can be extended straightfor- 
wardly to qudit (d-dimensional) systems where the indi- 
vidual detection error matrix D in Eq. (I) is replaced by 
a,dxd matrix. For independent detection of n-qudits, the 
overall error matrix M still has the tensor-product struc- 
ture as shown by Eq. (3), which allows easy calculation 
of M" 1 from D~ x . 

With the inverse error matrix M _1 , it is straightfor- 
ward to reconstruct the true distribution gi from the mea- 
sured data /j. The price we need to pay is that compared 
with A/j = \J fi(l — fi)/N, there is an increase of the 
standard deviation (error bar) Agi in our estimate of gi 
by the formula 9l = ^ 2 =1 M^ 1 fj. With some tedious 
but straightforward calculation, we find 



A 9i = [^(M^rfj-gfyN 



(8) 



As M _1 = ®fe = i -Dfc 1 an d 1 has matrix element 
1 — p' Q rs e p > 1 (when po ~ p\ ~ p <C 1), A'/ -1 has ma- 
trix element ~ e np which leads to exponential increase 



3 



of the error bar Agi with the qubit number n. To main- 
tain the same error bar Agi, the number of repetitions N 
of the experiment eventually needs to increase exponen- 
tially with n. For practical applications, this exponential 
increase of N by the factor e™ p is typically not a prob- 
lem for two reasons. First, as the detection error rate 
p is usually at a few percent level, the exponential fac- 
tor e np remains moderate even for hundreds of qubits. 
Second, this exponential increase only applies when we 
need to measure each element of the distribution <?j. In 
most of quantum information applications, we only need 
to measure certain operators which are expressed as ten- 
sor products of a constant number of Pauli operators for 
different qubits. In this case, N does not have the expo- 
nential increase as we show now. 

Suppose we need to measure an operator O, which is 
expressed as O = ®£ =1 o'fc' c , where <j l £ k is a component of 
the Pauli matrices when fj,k = 1, 2, 3 or the identity oper- 
ator when jttfc = 0. The number of the Pauli matrices n p 
in the tensor product expansion of O is called the support 
of O. To measure the operator O, we choose the measure- 
ment basis to be the eigenbasis of <rj£* for the fcth qubit . 
In this measurement basis, O is diagonal with the matrix 
element O = ®^ =1 diag(<T^ fc ), where diag(a l l l k ) — [1,1] 
for /ife = and diag(CT^ fc ) = [1,-1] for pk = 1,2,3. 
Under the distribution gi, the expectation value of O 
is given by (o) = ZA* = T.AT.jM^fj = 
Ej(Ei OiMf. 1 )^ = £\ 6p h where 6, denotes the di- 
agonal matrix element of O. Therefore, by defining a 
corrected operator O c , we can get the true expectation 
value (^O^ directly from the experimental data fj. Us- 
ing the relation M _1 = (S>fc=i DZ , O c is expressed as 
O c = <8>Li [diagHbW 1 ] • For Mfc = 1,2,3, 



diag(<7f )DZ 1 = [ 1 -1] 
= [ (1 - M) 

and for fik 



1-Po Pi 
Po 1-Pi 

-(1 - 2pi) ] 



(9) 



0, diag(a£ fc 



[1,1]. For simplic- 
ity of notation, we take po = pi — p. In this case, 
diag^)^ 1 = (1 - 2pr 1 diag( ( 7H for /i fc = 1,2,3, 
and the corrected operator O c is related with the origi- 
nal operator O by a simple scaling transformation O c = 
(1 — 2p) np O. The scaling transformation is indepen- 
dent of the qubit number n, so the error bar of docs 
not have exponential increase with n when the operator 
O has a constant support n p . 

The scaling transformation also applies to collective 
operators, but some caution needs to be taken for calcula- 
tion of their variance. For instance, if we take the collec- 
tive spin operator J z = Zk=i it is easy to see that 
J° = (1 — 2p)~ 1 J z as each of the terms of J z has support 



has non-uniform support for its superposition terms, 
one finds that (jf) c = n/4 + (1 - 2pf 2 (j 2 - n/4) = 
(1 — 2p) 2 \J 2 — (1 — p)l . With this transformation, 
we can correct the distortion to the spin squeezing pa- 
rameter by the detection error. Assume that the mean 
value of (J) is along the x-direction with (J) = (J x ) and 
the squeezing is along the z-direction. The squeezing pa- 
rameter is given by £ = \J n ( JJ) / ( J x ) 2 Using the 
transformation for (J 2 ) C and J^, we find that 



^n{{Jl) c )/{jtf = Je 



f 2 



(10) 



where — n 2 p(l—p)(l—2p)~ 2 (J x ) 2 is the contribution 
to £ 2 by the detection noise. After correction of the de- 
tection error, £ c gets significantly smaller compared with 
£ in particular when the qubit number n is large, and thus 
can be used to verify a much bigger entanglement depth 
using the criterion in Ref. @]. From Eq. (10), we find 
that the variation A£ c /A£ = £/£ c . As typically £ > £ c , 
the error bar for £ c after correction of the detection error 
gets significantly larger, and we need to correspondingly 
increase the rounds of the experiment TV to reduce the 
statistical error. 

To illustrate application of the error correction method 
here, as an example, we apply it to detection of genuine 
multi-partite entanglement in graph states. For a graph 
state \G n ) of n qubits associated with a g-colorable graph 
G, the genuine n-party entanglement can be detected 
with the following witness operator Q 



Wg„ =31-2 



e n (^+ i )/ 2 



1=1 \keQi 



(11) 



= 1. 



However, as J 2 = n/4 + J2k^i a k a i /4 which 



where Qi denotes the set of qubits with the Ith color 
(I = 1, 2, • • • , q), I is the identity operator, and Sk is the 
stabilizer operator for the A:th qubit (which is a tensor 
product of the Pauli operators for the feth qubit and 
(7^, for all it neighbors k' in the graph G). A state p has 
genuine n-qubit entanglement if tr (pWa n ) — (Wc B ) < 0. 
For an ideal graph state, all its stabilizer operators Sk 
have expect values (Sk) — 1- With detection error, the 
value of (Sk) gets significantly degraded. As an exam- 
ple, Fig. 1 shows the values of all (Sk) for two particular 
2-colorable graph states: a 10-qubit GHZ state (GHZio) 
and a linear cluster state (LCio). We assume 3% detec- 
tion error with po = p\ = p = 0.03 for each qubit. With a 
known magnitude p, the detection error can be easily cor- 
rected by a scaling transformation S^ = (1 — 2p)~ 7lpk Sk, 
where n p k is the support of the corresponding stabilizer 
operator Sk- Fig. 1 shows that after error correction, 
(S^) is almost unity. Its error bar increases a bit after 
the correction, but is still small. To show the influence on 
the entanglement detection, we assume the experimen- 
tally prepared graph state p ex corresponds to the ideal 
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Figure 1: Values of stabilizers before (lower points) and after 
(upper points) correction of the detection error (with the error 
rate p = 0.03) for the 10-qubit GHZ state (GHZio) and the 
linear cluster state (LCio). Error bars account for the statis- 
tical error by assuming N = 5000 independent measurements 
in each detection setting. 



Figure 2: The entanglement witness (Ws n ) under differ- 
ent state preparation errors p n for GHZ (GHZio) and clus- 
ter (LCio) states before (upper lines) and after (lower lines) 
corection of the detection error (error rate p = 0.03). The 
error bars are obtained by assuming N = 5000 rounds of 
measurements in each detection setting. 



target state pid distorted by small depolarization noise in- 
dependently acting on each qubit, so p ex = $ (pid), where 
the noise super-operator $ = <3)fc=i $fc an d Sfc {pid) — 

(l-3p n /4)p id +p„/4E ra)Z 4 P) ^4 /i) 0- InFi S- 2 - 
we show the witness {Wg„ ) as a function of the prepara- 
tion error rate p n , both before and after correction of the 
detection error (with an error rate p — 3%). For both 
GHZio and LCio states, without correction of the detec- 
tion error, we cannot detect any n-qubit entanglement 
even for a perfectly prepared state with p„ = 0. After 
correction of the detection error, we can confirm genuine 
n-qubit entanglement as long as the preparation error 
Pn i$ 5%. So, correction of the detection error signifi- 
cantly improves the experimental performance, and the 
improvement gets more dramatic when the qubit number 
increases. 

Finally, we briefly comment on the sensitivity of our 
error correction method to calibration of the detection 
error. In this method, the error magnitude p (or mag- 
nitudes pi, i = 0, 1, • • • , for general cases) is assumed 
to be known. If we have a relative error e in calibra- 
tion of the magnitude p, .i.e., Sp/p ~ e, the scaling 
transformation on the detected operator O leads to an 
relative error in the observed quantity 5 (^O^ / (o\ ~ 

2n p Sp(l — 2p) 1 ~ 2n p pe. As long as 2n p p ^ 1, which 
is typically the case as p <C 1, the relative error actu- 
ally gets reduced and the method here can tolerate some 
uncertainty in calibration of the error magnitude p. 

In summary, we have shown a method to correct any 
detection error through simple processing of the exper- 
imental data. The method applies to measurements in 
general many-particle settings, with or without separate 
addressing. Different from the conventional quantum er- 
ror correction, this method does not require encoding or 
change of the quantum circuits and works under arbitrary 
magnitudes of the detection noise, as long as the error 
magnitude has been calibrated. The cost of this method 



is moderate as it only requires repetition of the same 
experiment by some additional rounds to gain enough 
statistics and thus the method can readily apply to many 
experimental settings. 

Appendix:Proof of Eq (7). We can relate the L 
matrix to the M matrix defined in Eq. (3). Denote the 
space of n-bit binary strings with i bits of 1 as Si, and Si 
has dimension (™J . The matrix element M ap represents 
the probability of recording a n-bit binary string p as a, 
and Lij is the probability of recording a signal p G Sj as 
any string in the Si space. As a collective measurement 
does not distinguish the binary strings in the same space 
Si, Lij is related to M ap by L i3 = J2aeSi Mrp-The prob- 
ability Lij is apparently independent of the exact form 
of p, as long as p belongs to the space Sj, so we can 
pick up any p G Sj in Lij = ^2 aeS . M ap without alterna- 
tion to the result of summation. From Eq. (4) , we know 
M ilv = M fj.v (PqjPi). Let us define 

N jk ee £ A/~; = M pu (p' ,p[) = L jk (p' ,p[) , 

where v is an arbitrary element in S k ■ Now we show that 
N gives inverse of the matrix L: 

<?eSi j ueSj <?eSi 

In the second line, we have changed the subscript p in 
M ap to p as both /?, p belong to Sj. This proves Eq. (7) 
in the text. 
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